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Using Andreev and Lifshitz's supersolid hydrodynamics, we obtain the propagating longitudinal 
modes at non-zero applied pressure P a (necessary for solid 4 He), and their generation efficiencies 
by heaters and transducers. For small P a , a solid develops an internal pressure P ~ P%. This 
theory has stress contributions both from the lattice and an internal pressure P. Because both 
types of stress are included, the normal mode analysis differs from previous works. Not surprisingly, 
transducers are significantly more efficient at producing elastic waves and heaters are significantly 
more efficient at producing fourth sound waves. We take the system to be isotropic, which should 
apply to systems that are glassy or consist of many crystallites; the results should also apply, at 
least qualitatively, to single-crystal hep 4 He. 
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I. INTRODUCTION 

In 1969 Andreev and Lifshitz developed a theory of 
supersolids. 1 Although the microscopic physical descrip- 
tion was for flow of vacancies, the macroscopic equations 
did not depend on vacancies in an essential fashion. At 
about the same time Thouless 2 - and Chester— both sug- 
gested the possibility of superflow in a solid by vacan- 
cies. In addition, Leggett4> pointed out the possibility 
of Non-Classical Rotational Inertia (NCRI) associated 
with quantum-mechanical flow via a superfluid velocity 
(a phase gradient) opposite the local velocity of the ro- 
tating lattice. 

Since the observation of NCRI by Kim and Chan^i 6 - a 
number of laboratories have reproduced their worki 7 - - — 
(For reviews that emphasize experiment, see Refs.lTsI and 
Ho .) Were NCRI the sole criterion for superflow of solids, 
there would be strong reason to accept that such super- 
flow has been observed. However, a supersolid should 
also have other properties, including a fourth sound-like 
mode, as predicted by Andreev and Lifshitz, and modi- 
fled elastic waves with higher velocities, since the super- 
fluid mass does not participate in the motion. (We re- 
mind the reader that a fourth sound mode in superfluid 
He occurs only when the normal fluid is entrained by a 
porous medium; in the present case the lattice serves as 
the porous medium.) Neither a fourth sound mode nor 
velocity shifts have been observed.— ~— However, a stiff- 
ened shear response is observed ; 20 ' 21 although not enough 
to explain the observed NCRI. 16 Note also recent work 
indicating that supersolidity in 4 He can only occur below 
55 mK.^ 

As a guide to experiments to observe the fourth sound 
mode^ the present work calculates various quantities 
relevant to its observation, such as the relative efficien- 
cies of a transducer and a heater in producing both lon- 
gitudinal elastic waves and fourth sound waves. It also 
considers the effect of a non-zero applied pressure P a ; to 
solidify 4 He, even near T = 0, requires P a > 25 bar. 
To our knowledge, previous works have not included the 
effect of P a . 



Although we believe that vacancies are essential to 
a microscopic understanding of superflow in solids, in 
the hydrodynamic theory they play no fundamental role, 
other than as an additional variable largely tied to diffu- 
sion. Indeed, we believe that the hydrodynamic theory is 
more likely to describe a supersolid related to the NCRI 
effect proposed by Leggett than to vacancy superflow. 

Most of the present work assumes that the system is 
isotropic. One effect this has is that the superfluid den- 
sity, which properly is a second rank tensor p^ , is propor- 
tional to the unit matrix, so we take p s « 1 p 8 , 23 i 24 We 
then write the superfluid fraction as 



where p s is the superfluid density and p is the total (mass) 
density. f s is unity in a superfluid at low temperatures. 
However, in putative supersolid 4 He, the measured NCRI 
fraction, which if due to superflow should be equated to 
f s , is never greater than about 0.2. The effective nor- 
mal fraction /„ thus has the curious property of being 
not less than 0.8, although at T = there are no excita- 
tions to destroy the superflow. We have previously noted 
this difficulty, 25 and proposed that the lattice be given a 
mass fraction in addition to a contribution f^ x due 
to excitations, so that 1 = f s + f^ x + ft,. This permits, 
at T = 0, no excitations (so f% x = 0) but f s < 1. In 
this viewpoint, the lattice velocity is identified with in, 
where m is the lattice displacement, and the effective 
normal fluid fraction /„ is the sum of fr, and a part ff x 
due to excitations: /„ = f% x + fi,. 

It is known that the more annealed (and thus more 
crystalline) the sample of 4 He, the smaller the NCRI 
fraction.! 1 Likewise it is known that the more quenched 
the sample, the larger the NCRI fraction* 1 ^ Hence the 
supersolidity is more likely to occur for less crystalline 
samples, which might be either glassy 2 ^ or consist of a 
large number of small crystallites.— In both of these cases 
an acoustic probe is likely to take a rotational average, 
thus making the system behave more like an isotropic sys- 
tem than a crystal. Therefore we consider systems whose 
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macroscopic properties are isotropic. If a pure crystal of 
hep 4 He were to be supersolid, then the results we obtain 
will be only an approximation; nevertheless they will be 
a useful guide for experiment. 

We also note that we are working in the linear regime, 
where the disturbances produced by a heater or trans- 
ducer are expected to be only a small perturbation, as 
is assumed in all theories of this sort. It is possible that 
solid 4 He is ultrasensitive to temperature or to stress (e.g. 
if it is a glass, perhaps the atoms can be driven off their 
sites by a transducer). Nevertheless, there should always 
be a linear regime; Ref. [l9| notes that their membrane for 
producing putative fourth sound produced strains much 
below the critical value. 

The present work is intended to be self-contained, 
although we do not derive the equations of motion, 
which may be most explicitly obtained from Sect. IV of 
Ref. I28li22 nor do we derive the P a dependence of various 
quantities, which is done explicitly in Ref. l30l . 

Sections |TT] and [TTT] present the thermodynamics and 
equations of motion, respectively. Section HVl studies the 
eigenfrequency and eigenmodes for each of the longitudi- 
nal propagating modes. Section IVl finds, for each mode, 
the stress and temperature response in terms of the nor- 
mal fluid velocity. Section [VTI finds and discusses the ef- 
ficiency of generating each propagating mode by a trans- 
ducer and a heater. Section TVIII provides a brief summary 
of the results. 

One of the features of the hydrodynamic theory of su- 
persolids is that it contains stress due to both internal 
pressure P and the lattice, in order to permit the system 
properties to continuously transform into those of a su- 
perfluid. This means that under an applied pressure P a 
the pressure and the lattice each take up a part of it. We 
believe that use of both a pressure and a lattice stress is 
needed not merely for solid 4 He but for other solids as 
well, particularly those under pressure or with point de- 
fects that are not in equilibrium^ Appendix [Al discusses 
the relationship between the internal pressure P and the 
applied pressure P a , which we estimate using experimen- 
tal data and results from Ref. [3(| Appendix [B] discusses 
the relative sizes of velocities and strains in a crystal un- 
der an applied pressure. Appendix [D] finds the relative 
size of two thermodynamic derivatives of temperature T 
that appear in the generation efficiencies. 



II. THERMODYNAMICS 

The thermodynamic equations for a supersolid are 
given in terms of the energy density e, entropy density s, 
unsymmetrized strain = diUj , mass density p, super- 
fluid velocity v s , and momentum density 

g = p n v n + PsV Sl (2) 



and their thermodynamically conjugate quantities 
Here v n is the normal fluid velocity. Specifically, 

de = Tds + \ lk dw ik + pdp + v n ■ dg + j s ■ dv s , (3) 
e = -P + Ts + XikWik + PP + v„ ■ g + j s ■ v s , (4) 
= — dP + sdT + WikdXik + pdp + g ■ dv n + v s ■ dj 8 . 

(5) 

Here the thermodynamically conjugate quantities are 
temperature T, (unsymmetrized) elastic tensor density 
Xik (with units of pressure P), chemical potential p (with 
units of velocity squared), normal fluid velocity v n , and 

j s = g- pv n = p s (v s - v n ). (6) 

Note that j s has units of momentum density but is in- 
variant under Galilean boosts. That is, if both v n and v s 
are boosted by Sv, j s does not change. 
We find it convenient to define 

j n =pv n , (7) 

so that 

g = jn+3s- (8) 

Unlike j s , the quantity j n is a momentum density both in 
units and in its properties under Galilean boosts; under 
a boost by 5v, both g and j n are boosted by p(Sv). 

III. HYDRODYNAMIC EQUATIONS 

Eq. ([3]) shows that there are five independent thermo- 
dynamic variables. Two of them are scalars (s and p), 
one is a tensor (wtk — diUk) and two are vectors (g and 
v s ). In developing the hydrodynamic equations we will 
employ the first three, but we will use the two vectors j n 
and j s in place of g and v s . For an ordinary solid, where 
v s does not appear, it is convenient to use the variables 
a = s/p, p, Wik, and g, since a decouples from the other 
variables. Such decoupling does not occur for the super- 
solid. Note that one could also use the scalar variables T 
and p, or T and p, or s and p. 

Unless otherwise specified, thermodynamic derivatives 
with respect to p, s, or Wik are taken with the other two 
variables held constant. 

We consider small amplitude excitations of the form 
exp [i(k ■ r — u>t)], where the wavevector k is taken to 
be known. Then, with primes denoting deviations from 
equilibrium, in the absence of damping the equations of 
motion are given by^i^ 



p' + d % g[ = 0, (9) 

g'i - d k a' ik = 0, a lk = Xik - PSik, (10) 

v' si + dip' = 0, (11) 

s' + Sd^ = 0, (12) 

< - Vn', = 0. (13) 
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Rather than the stress tensor the momentum flux 
n ife = -a ik + p s v si v sk + p n v ni v nk ~a ik has also been 



1.31.32 



employed,— ^ as well as gi = j, 

In terms of the thermodynamic variables p and s, we 
have 



p 



dp , df£ s , dp w , 
dp ds dwji J '' 



da 



ik i 
S 



d(Tik 



(14) 
(15) 



lK dp ^ ds dwji jr 
The equations of motion (jSJ), (|T2"j). and (fT3")) directly give 



P = 



hi Vm. a 



(16) 
(17) 
(18) 



The other two equations of motion, (fT0| and (fTTj). can 
now be written in terms of v' si and u^, or, equivalently, 
j'si and 

From (fTK)) - (fT8|) . momentum conservation (fT0|) gives 



kkkjp~ 1 (daik/dwji) in (|20l) contains two terms, one pro- 
portional to <5ij and one proportional to kiki. 

Taking the dot product of these equations with fcj then 
gives two equations in the unknowns u> 2 , k-j' s , and k- j' n . 
This yields w 2 and the ratio k-j' s /k-j' n . In addition, taking 
the cross-product of (j2"2"j) with fcj gives, since fc x u s = 0, 
identically zero. Further, taking the cross-product of ([20)1 
with fci gives an equation having terms proportional to 
ui 2 and k 2 , both multiplying k x v n . 

There are two ways to solve the resulting equations for 
k-j' s , k-j' n , and kxv n . One solution is to take k x v n = 
and kij' sl ^ 0, kij' nl ^ (purely longitudinal modes), 
with the frequencies determined by the two equations 
in the unknowns cj 2 , k ■ j' s , and k ■ j' n . Hence this set 
of modes is purely longitudinal. The other solution is to 
take kxv n ^ and k-j' s = 0, k-j' n — (purely transverse 
modes), with the frequencies determined by the the cross- 
product of (|20p with fcj. v s does not participate in the 
transverse modes, so their mass weighting involves only 
p/pn, and their frequencies squared should be higher than 
in the normal solid by p/p n - To our knowledge such an 
effect has not been observed^ 



ufA = - k k a' tk 

_ k k ki da ik 
ui dp 



■9i 



kkk da ik 

Ld ds 



k k kj d(7i k , 
a Vn i- 



(19) 



We now rearrange to use the variables j' n and j' s . From 
©-©, multiplying (JJ]) by uj gives 



IV. LONGITUDINAL EIGENFREQUENCIES 
AND EIGENMODES 



Recall that, unless otherwise specified, thermodynamic 
derivatives with respect to p, s, or Wi k are taken with the 
other two variables held constant. 







k k ki- 



do;,, 



dp 



Jsl 



L0 2 8 a + ( fcfefc, ^ 



-\-k k k\ 



s daj k 
p ds 



k k kj 



a -t- \ ki, 

1 dajk 
p dwji 



dp 
Jn'f (20) 



Likewise, from (1141) and (1161) - (1181) . the superfluid equa- 
tion of motion (fTTj) gives 



kiki dp , kiki dp 



ui dp 



fli 



uj ds 



Ml 



U! dlVjl 



■Jnl- 



From ©-(El) and p s v s [ = j, 
(PHI) by p s uj and rearranging gives 



(21) 

(Ps/p)jn'i, multiplying 



js'i + fs 



uj 2 5 z [ 



u 2 5 a - kihf s p^ 

— p (kiki— — h kiki — kikj — — 

\ dp p ds p dwji 



jn'f (22) 



Equations (|2"0|) and (|2^|) yield the normal mode fre- 
quencies and their eigenvectors (the ratio of the responses 
of the normal and superfluid currents). In what follows 
we consider only an isotropic solid. The effect this con- 
straint has on (|2"U)l and (j2"2"j) is that the second-rank ten- 
sors are all proportional to the unit tensor, and the term 



A. Some Properties and Definitions 

We now compute the quantity dai k /dp, which appears 
in (|2H)) . We take the strain response of a solid to P a to 
be isotropic (i.e., ~ 5i k w^\ where the superscript 
(0) denotes the static value). Recall that Wi k is unsym- 
metrized; here we take only the static part, due to P a , to 
be symmetric, as does Ref. [34]. Then, by Ref. H3, 



A 



(0) 
ik 



K 



\ x C°) i ( (°) i (°A 

Pv I 5 lk w\ l + p v (w 2 y + w ki 'j 



so that we can write 



d\ r 



OX 



dp dp 



x (°) 
(23) 



(24) 



Here, K and pv are the respective bulk and shear mod- 
uli, with units of P; py is completely distinct from p. 
Eq. (|24"|) is also employed in Ref. [30, although there 
a = s/p is held constant rather than s. At T ss 0, the 
difference should be negligible. Thus we can write 



da 



ik 



d\ik 



dp dp 



Si, 



dP 
dp 



dX dP 
dp dp 



da 

~d~ P 



Sik, 



(25) 
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where we use a to distinguish a stress (with the same 
units as a^) from a = s/p; a and a are not related. Note 
that da /dp is not a true derivative, merely a definition; 
further, we do not here define a a. 

We now compute the quantity daik/dwji, which also 
appears in ([20]) . Since Ref. |30 shows that (dP/dwik) ~ 



,(0) 



v% 5ik, we can write 



dP _ dP 

OWik OW 



(26) 



dP/dw is evaluated in Ref. yO, and is given in Ap- 
pendix [B] We also use the definitions 



OX T ^ 4 



9a_dX_dP 

dw dw dw ' 



As above, da /dw is not a true derivative, merely a defi- 
nition. Eqs. (|23j and ([26]) then give 



~d~uy = (^w _ ■ " /,S '''''''' ''' + /M 



(28) 



Thus, 



da ik ., 



92F 



Wi^f/m = [q^- Hk ■ 3n) + Vvk 2 f ni . (29) 
For fc • j' n ^ and k x j„ = (the longitudinal case), 

U2 ■/ 



(30) 



h(k ■j' n ) = k 2 j' ni so that JHD gives 



It is convenient to define the velocities cq and c\ , which 
satisfy 



2 _ dp 

Co=/ V 



da 1 9<T 
<9p p dw 



(31) 
(32) 



If (j, rather than s, were held constant, then Co would 
be the sound velocity in an ordinary fluid, and c\ would 
be the velocity of sound in an ordinary solid with no 
supcrflow. 30 Using the Gibbs-Duhem relation ([5]) and ne- 
glecting thermal expansion and terms second order in ve- 
locities gives 



8X 



dP 2 

9^ = Co + ^ dp 



Then, eq. (|25j) gives 

da 
dp"* 



(i-4 0) ) 



dX 



(33) 



(34) 



In the following we use either (|25l) or (|3"4"|) . depending on 
convenience. 



B. Reducing the Equations of Motion 

Momentum Equation ([20]) : We take s(daik/ds) —> 
0, which should be a reasonable approximation for solid 
4 He at low temperatures, both because s — > as T — > 0, 
and because K and py (and therefore Xik at constant 
Wji) should be nearly independent of s. Substituting (|25[) 
and (30|) into ([20]) and using ((32]) and ((34]) then gives, for 
a purely longitudinal mode, 







~2 , (o) dX . ., 
C + 11 dp ' 



■ ' i r 2 2 7 21 • / 

J s +[U - C x fc J J n ; 



(35) 



~2 _ 2 
C = C n 



(36) 



here we define, to simplify the equations, 

dX 
dp 

Appendix [B] finds that \dX/dp\ ^> c 2 ., so that c 2 « 
—dX/dp. It also finds that c 2 is expected to be posi- 
tive, and first order in P a /K. Further, it shows that for 
P a < K we have c\ > c 2 » eg. 

Super fluid Equation (|2"2"1) : A Maxwell relation that 
follows from ([3]), combined with ([Ml) , gives 



9w 



. dX 



<9p 9p 



(37) 



Then, neglecting s(dp/ds) = s{dT/dp) ~ T 4 , and taking 
the mode to be purely longitudinal, eq. ([22]) gives 



= [c^ 2 - / s c 2 fc 2 ] + /, [w 2 - fc 2 - 



(38) 



We use ([35l and ([38]) first to find the longitudinal mode 
frequencies, then to find the superfluid-to-normal ratios 
of current density and velocity in each longitudinal mode. 
For f s -> 0, eq. fl38) gives either w 2 = fsC^k 2 (fourth 
sound) or j' s = (no superflow). In the latter case, sub- 
stitution into (j3"51 then gives oj 2 = c\k 2 (first sound). 



C. Longitudinal Mode Frequencies 
Eqs. ([35]) and ([38]) yield 
Q=lo\1- f s )-u 2 k 2 



c\ + / s c 2 - f s ( 2c 2 + u>„ 



:0)dX 

dp 



2 2 ^2 / ^2 

C 1 Ln 



c + w 



(0) 



Solving ([39]) to first order in f s gives 



(39) 



k 2 ~ k 2 ~ /s 



(40) 



and 



^,2 - y. - J* 1 c o 



,(°) 



72 



9Ac 
9p c 



/? ~2 



(41) 
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In the limit where Cg <C dX/dp and w 
K), 



(0) 



<C 1 (i.e., P a « 



(42) 



With w2 ~ for /, <£ 1 



Js4 
in4 



c l 



Max 

V U dp 



(48) 



Appendix iBl finds that both terms on the right-hand-side 
of (|42p are second order in P a /K. Further, it shows that 
for P a -C K we have c 2 3> c 2 » c 2 ,. 



(o) 



Appendix IBl shows that if P a /K <C 1, then w 
and c 2 , <C dX/dp. Thus, eq. ([56"]) gives c 2 « —dX/dp > 



< 1 



,(o) 



dA/dp. Then, 



D. Longitudinal Mode Structure — Currents and 
Velocities 



We now find the ratios of the normal fluid and super- 
fluid response for both longitudinal modes. These ratios 
will be used to calculate, for each mode, the response to 
the stress and temperature produced by transducers and 
by heaters. We employ 



P PsW s -v' n ) 



Pa 



P< 



1 = 14 

fs Jn 



1. 



(43) 



The ratios j' s / j' n for each mode can in principle be ob- 
tained from the normal mode frequencies and either of 
©5} or (US}. 



1. First Sound Mode Structure 
From (|38|). with the subscript 1 denoting first sound, 

O'.l 



_ f k 2 '- 



(44) 



JS4: 

jn4 



c l 



(49) 



Then, using (1431) . the ratio of superfluid velocity to nor- 
mal velocity for fourth sound is, to lowest order in / s , 



u s4 



fs? 



(50) 



Appendix [B] shows that c\ 3> c 2 , 



so > w„ 4 . 



V. LONGITUDINAL MODES - STRESS AND 
TEMPERATURE RESPONSES 

We now calculate the deviations from equilibrium of 
the longitudinal stress and temperature produced by a 
transducer and by a heater. We consider that only the 
a'n component of the stress is generated. For notational 
simplicity we employ &' = o' Y \\ recall that a is reserved 
for the entropy/mass. 



Stress 



Substituting w 2 from (j40|) . accurate to zeroth order in / s , 
gives a ratio accurate to first order in f s : 



4 1 

Jnl 



-fs 1 



(45) 



Then, using (|43|) - the ratio of superfluid velocity to nor- 
mal velocity for first sound is 



Conservation of momentum (|T0|) yields 



Substituting the ratio j' s /j' n from (|45|) and (|49|) and 
Wl 4 = +c 1A k from (gDJ) and flU} into ([ST} gives the 
stress associated with each mode. For / S <1, 



si 



Appendix [B] shows that c 2 » c 2 , sov^ » . 



(46) 



Cl 



0-4 ~fs ^pC V ni 



Pv'nl 



fi^pcQ—v. 



-pciv nl , 



nV 



(52) 
(53) 



2. Fourth Sound Mode Structure 



where we have used c 2 3> c 2 (see Appendix [B]) . 
The total stress deviation therefore is 



From (|35|). with subscript 4 denoting fourth sound, 



k 2 



3ni 



Max 

"ll dp 



(47) 



-cipv nl 



1 - tfco 



Ci V, 



(54) 



The ratio v' n4 /v' nl depends on the mode generator, to be 
discussed in the next section. 
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B. Temperature 



A. Transducer Properties 



The temperature deviation is less straightforward to 
obtain because it is a function of the variables s, p, and 
w 3l : 

dT , 



dT 8T , 
os dp 



dw 



(55) 



Since dT/dwji = dXji/ds, and K and pv depend only 
weakly on s, by ([2"3"]) we neglect dT/dwji. Substitution 
for p 1 and s' from (fT6|) and (| 1 T|) then yields 



dTki ,dTh , 
T ps s- u ni + g v 

OS W Op L) 



(56) 



We earlier showed that the mode is longitudinal, so we 
drop the indices i. The identity g' — p (1 + j' s / j' n ) v' n 
then yields, for both modes, that 



r 



dT dT 



ds 



P 



dp 



i 



J'n 



(57) 



Substituting the ratio j'Jj' n from (g5J and (gSJ) and 
wi 4 = +cfci 4 from (l40l) and (|4T|) into ([57]) gives the tem- 
perature associated with each mode. To lowest order in 
f s we obtain 



■'n 1 



Cl 

Co 



<9T <9T 

S— h p 



ds dp 
dT 
ds 



dT 
dp 



1 



(58) 



(59) 



For a solid at low temperature, Appendix [C] gives 
pidTjdp) ps lQs{dT/ds). In_addition, for P a < AT (as 
is the case here), Appendix 
(ESI- (EH become 



gives c 2 ;§> c 2 . Therefore 



T i ~P^{ 1 

dp a 



n 



_ f -\ST±vA__ f - 
dp c £ cq 



dr c| v n ' 4 v n [ 

1 dp c 2 v n [ c 



The total temperature deviation therefore is 



T' = T[ 



n 



dp ci 



i - /; 



C^CO Vn'i 



(60) 
(61) 

(62) 



The ratio v' ni /v' nl depends on the mode generator, to be 
discussed in the next section. 



VI. LONGITUDINAL MODE GENERATION 

A transducer produces, and therefore can be used to 
detect, stress deviations. A heater produces, and there- 
fore can be used to detect, temperature deviations (when 
used as a detector, a heater is called a thermometer). To 
utilize the results of Sec. [V] we find v' nA /v' nl for each 
device, then substitute it into ([52 ]) -([51 ]) and ([50 ]) - ([52" ]) to 
find the respective stress and temperature deviations pro- 
duced by transducers and heaters. 



For a transducer we take v' s = v' n (and therefore j' s = 0) 
so that 



Q = j' s \txn = Kl+isJtrn 



hi j' -I- h±n' 

■I Jnl~ r •/ Jni 
Jnl JnA 



(63) 



where the subscript "trn" denotes properties of a trans- 
ducer. Then 



Jni 
n Jnl 

Use of gSJ and (gSJ) yields 



Js l/ Jnl 
JsaI JnA 



u n4 



C 2 / C 2 



fs? 



(64) 



(65) 



where we have taken cf ^> c 2 (see Appendix IB"]) . 

Eq. ([52} gives in terms of u^, regardless of gener- 
ator. Use of (JB5J) in ([5"3]) - (j53]l gives 



4 I trn 



fs 2 pC v' nl , 



pciv nl 1 + /.. 



Cl 



-pciu nl , 



(66) 
(67) 



for /„ C 1 and c 2 , <C c 2 . Thus the stress produced by 
a transducer primarily goes into first sound, with a frac- 
tion f^ 2 (co/ci) of the stress going into fourth sound. 
Eqs. ([52"]) and ([55|) divided by ([5T]) are the two entries in 
the top left of Table [H 

Eq. ([60]) gives T[ in terms of v' nl , regardless of gener- 
ator. Substituting (155]) into (15TT) yields 



T'l 

J 4ltr 



J Js P „ — ■ 

dp c 



(68) 



By flU]), T4 = fl /2 ( Cl /co)T{, so T x ' and T4 could be of 
the same order of magnitude. Eqs. ([60]) and ([68]) divided 
by ([67]) are the two entries in the top right of Table Q] 



B. Heater Properties 

For a heater we take g' = 0, so that 

=SV - [(jii + Jnl) + C?i 4 + Jn 4 )]htr 

'j. 



^ + 1 ) Jnl I htr + ( '^T + 1 ) in4lhn • 
Jnl / \Jn4 



(69) 



where the subscript "htr" denotes properties of a heater. 
Then 



J n4 



htr 



Jra4 



htr 



(jsl/jnl) + l 

(jVi;j+i' 



(70) 
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and substitution from (1451) and (|49[) yields, for /,<!, 



u n4 



(71) 



htr 



Here, we have used cf 3> c 2 (see Appendix [Bjl . 

Eq. ([rJtJ|) gives T{ in terms of v' nl , regardless of gener- 
ator. Substitution of dTJ) into ([6"T ]) -([6"2" ]) gives 



J 4lhtr 



dp c 



hlr 



dp ci 



Co 



-1 



nl 



9/9 c 



(72) 
(73) 



for /, <C 1 and c§ <C cf. Thus the temperature pro- 
duced by a heater primarily goes into fourth sound, with 
a fraction f s (co/ci) of the temperature going into first 
sound. Eqs. dHOj) and ([72]) divided by <[73]) are the two 
entries in the bottom right of Table [I] 

Eq. (|52p gives cr^ in terms of v' nl , regardless of gener- 
ator. Substituting (fTTj) into ([53"|) yields 



Mlhtr 



(74) 



Since / s <§; 1 and c^ <C c 2 , we have ^ 0-4. Eqs. ([5J 
and (|74p divided by (f75|i are the two entries in the bottom 
left of Table H 



Although a transducer primarily produces stress, 
it also produces a small temperature deviation T'. 
M22/M12 gives 



Ti 



1/2 



£1 
Co' 



(transducer) 



(76) 



For / s < 1 and C 2 , <C c\ , it is not clear which of the terms 
in (|T6"|) dominates. 

A heater generates mostly temperature, T', and .M42 
shows that a heater is efficient as a fourth sound genera- 
tor. M.A2/M.Z2 gives 



/7 1/2 ^»i, 

Co 



(heater) 



(77) 



which is large. In fact, for f s ^l and c 2 , <C c\ , .M42 = 
T^/T' ?s 1 and therefore nearly all of the temperature 
response corresponds to the fourth sound mode. 

Although a heater primarily produces temperature, it 
also produces a small stress deviation 
gives 



1/2 £0 
ci ' 



(heater) 



(78) 



Therefore, with f s <C 1 and eg < c(, eqs. ([75j and ([78)) 
imply that stress deviations do not contribute an appre- 
ciable amount of fourth sound, whether produced by a 
transducer or a heater. 



Generation Efficiencies 



VII. SUMMARY 



A proper treatment of the response of a given detec- 
tor (transducer or thermometer) to a given mode (first 
or fourth sound) would consider the incoming mode 
and what happens under reflection from the detector; 
this would give the net stress and temperature at the 
detector*^ We consider only the issue of generation. 

A transducer generates mostly stress, a' . Thinking of 
the equation entries in Table [J as a 4-by-2 matrix Ai, 
M.\\ shows that a transducer is efficient as a first sound 
generator. M-ixj M.\\ gives 



We have studied the implications of the Andreev and 
Lifshitz theory of supersolids for the generation of a 
fourth sound mode in a solid under an applied pressure 
P a , including the relative efficiencies of a transducer and 
a heater in producing both longitudinal elastic waves and 
fourth sound waves. The present results apply when the 
bulk modulus K 3> P a . 
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Heater 


1st Sound 


7fJ~fs 2 ClC — 


/s CO 


4th Sound 


rpi ~ ,/s C — Js ^ rp f <K. 


T| 1 

Cl 
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Appendix A: Relating Applied Pressure P a and 
Internal Pressure P 



We now use the experimental data of Ref. |35j to esti- 
mate P/P a - We then evaluate when P a P. 
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Ref. [30j gives that 



P 

K 



K P a 

2K 2 ' 



where A* = A - V(dK/dV) 

k (0) 



Wik ,f,7V ' 



(Al) 

Under hydrostatic 



compression X\J - P8 l k = -P a §ik- Thus, 



A[ ; 0) = 3P - 3P a 



3A ^ 
2A 2 



-3P«. 



(A2) 



Unfortunately, A* is not a quantity measured exper- 
imentally, since the structure of the energy density dic- 
tates that the derivative is taken at constant strain (i.e., 
constant lattice site density). 

In what follows, we roughly estimate dK/dV by as- 
suming it to be of the same order of magnitude whether 
taken at constant Wn- or under typical experimental con- 
ditions. That is, we take 



dK 
dV 



dK 
dV 



cr,7V,exp 



AK 
~A~V 



(A3) 



(T, 7V,exp 



We now consider the data of Ref. [35J. Although the 
samples were necessarily under pressure, Ref. HH appears 
to apply c\ = (l/p)(dX/dw) = (l/p)[K + (4/3)//y] with- 
out including corrections due to P a . 30 Nevertheless their 
result should permit a rough estimate (for simplicity we 
consider that T sa 0). We use K = (1/3) [cu + 2c 13 ], 
where cu and C13 are elastic constants^. Select parts of 
Tables I and II of Ref. [35J are reproduced in Table [TT] for 
two molar volumes, which is sufficient to make estimates. 

These data give AV « —1.2 cm 3 /mole and AK « 
160 bars, so that AK/ AV w —133 bars- mole/cm 3 , which 
we take to be constant since the elastic constants in Fig- 
ure I of Ref. [35| are linear in volume. Thus we obtain the 
two values for K* in Tabled K* w 10K at P a = 31 bars 
and K* 6.6A at P a = 52 bars. 

For P a « 31 bars, P/P a « 0.53 (and thus X n /P a « 
0.47). For P a w 52 bars, P/P a « 0.37 (and thus 
A n /P a « 0.63). 

For U = 20.5 cm 3 /mole, P a /A" « 0.10, so it is ap- 
propriate to take P a <C A. For V = 19.28 cm 3 /mole, 
P a /K r; 0.11, and P a <C A is still a reasonable approx- 
imation. Therefore, for applied pressures less than 100 
bars (and possibly higher), P a <C A likely holds. 

P a /K increases as P a increases. Therefore, although at 
higher P a we may find that P 3> P (extrapolating from 
Table [TTJ) , P a /K might become on the order of unity, 
and the approximations made in the present work and in 
Ref. no longer apply. 



Appendix B: Velocities, Thermodynamic Derivatives 
and Strain in a Crystal Under Applied Pressure 

We now estimate the relative sizes of the velocities c\ , 
Co, c, and Cq in the limit P a <§C A. We use the relation- 
ships between thermodynamic derivatives and applied 
pressure given in Ref. |30j. 



To lowest order in P a /K, Ref. [30 gives 



w u 


Pa 
A' 




d\ 


_ VP a dK 




op 


pK dV 


(7,W 


dP 


V 2 P 2 a d 2 K 


dp 


2pK 2 dV 2 


dP 

dw 


= -Pa (l - 


V 
A 




K+lp v 





p 



(Bl) 
(B2) 

(B3) 
(B4) 
(B5) 



where a — s/p. Here the internal pressure P has been 
taken to depend only on the square of the strain. Al- 
though dX/dp and dP/dp in Ref. |30j are taken at con- 
stant ex, not s, at solid 4 He temperatures we assume that 
a w w s. Note that to lowest order in P a /K, the strain 
of eq. (|B1[) agrees with Ref. H3, which includes lattice 
stress but no internal pressure P. (Ref. [3(3 also finds a 
P 2 term in the strain that is not obtained in Ref. l34l ) 
We now use the Gibbs-Duhem relation ([5]) to determine 



Cq. 



dp 
dp 



dP 
%~ 

' pK 2 



,(0) 



OX 



" dp 
V d 2 K 
~2 dV 2 



dK 
dV 



(B6) 



Note that Ref. Stakes p(dp/dp) = dP/dp, and thus does 
not include the term proportional to the static strain. As 
for P, p depends only on the square of the strain, via a 
Maxwell relation. This is not true for a good liquid. 
Eq. (|B6|) shows that c§ is second order in P a /K, 



whereas (|B2[) shows that dX/dp is first order in P a /K. 
Thus, for P a < A, \dX/dp\ > | eg | . Therefore §§§ gives 



? 2 



-dX/dp > c 2 . 



(B7) 



Further, we may find the sign of c 2 . A is a measure of 
the stiffness of a solid. Thus, as V increases at constant 
particle number and strain (or lattice site number den- 
sity), i.e., as vacancies and lattice sites are added to the 
system, A should decrease, or dK/dV < 0. Then (|B2|) 
gives dX/dp < 0, so that 



-dX/dp > 0. 



(B8) 



Moreover, since to first order c 2 is independent of P a , 
cqs. ||B2j) and ||B5j) give 



(B9) 
(B10) 



c\ > \dX/dp\. 
Combining JB7| and (JB9J) yields 



C? > C 2 > C„. 
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Volume 


Pa 


C33 


C13 


K 


AK/AV 


K* 


P 


P/Pa 


Pa/K 


(cm 3 /mole) 


(bars) 


(bars) 


(bars) 


(bars) 


(bars mole/cm 3 ) 


(bars) 


(bars) 






19.28 


51.6 T 


980 


198 


460 


-133 


3020 


19.1 


0.37 


0.11 


20.5 


31.4 f 


630 


142 


300 


-133 


3030 


16.6 


0.53 


0.10 



TABLE II. The first four columns are experimental data on hep 4 He from Ref . |35|. We also employ K = [C33 + 2ci3]/3, assume 
K to be linear in V, take K* w K - V(AK/AV), and find P/P a from $EQ. f Value was estimated by extrapolation from Table 
I of Ref. HI 



Finally, (j42|) gives c 2 , ~ c 2 — (c 4 /c 2 ), which implies that 
S 2 < S 2 . Thus, (IBTOI) gives 



c? > c 2 > c^. 



(Bll) 



The quantities Cq and Cq may be of the same order. 



Appendix C: On s(dT/ds) and p{dT/dp) 



For an insulating solid at low temperatures, it is well- 
known that s = ^(T/u) 3 , where £ is a dimensionless 
constant and u is the mean velocity of longitudinal and 



transverse ordinary sounds It immediately follows that 
dT 



T 
3"' 



and that (dT/du) 

dT 
^ dp 



T/u. We then have 



u dp 



dT 


du 


s ^ du 


s dp 



(CI) 



(C2) 



We are at low enough temperatures that we may consider 
u to depend only on density, so that both s and T may 
be considered nearly zero. 

We now estimate {p/u)(du/ dp)\r = —(V/u)(du/dV)\T 
of (|C2[) , where the molar volume V ~ p" 1 . Us- 
ing data from Fig. 13 of Ref. we take an averaged 
longitudinal sound velocity of 5 x 10 4 cm/s to be u, 
V to be 20 cm 3 /mole, and (from Fig.13) Au/AV « 
(du/dV) to be -0.83 x 10 4 mole/(cm 2 -s). This gives 
(p/u)(du/dp)\T ~ 3.3. Comparison with (jCip yields 
pidTjdp) « 10 sidTjds). 



